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Abstract 

We discuss central exclusive production (CEP) of the tensor Xc(2^) meson in proton- (anti) proton 
collisions at Tevatron, RHIC and LHC energies. The amplitude for the process is derived within 
the /ct-factorisation approach. Differential and total cross sections are calculated for several unin- 
tegrated gluon distributions (UGDFs). We compare exclusive production of all charmonium states 
Xc(O^), Xc{'^~^) and Xc{'^~^)- Equally good description of the recent Tevatron data is achieved both 
with Martin-Ryskin phenomenological UGDF and UGDF based on unified BFKL-DGLAP ap- 
proach. Unlike for Higgs production, the main contribution to the diffractive amplitude of heavy 
quarkonia comes from nonperturbative region of gluon transverse momenta Q± < 1 GeV. At y ~ 
0, depending on UGDF we predict the contribution of Xc(l^i2+) to the J/^ + 7 channel to be 
comparable or larger than that of the Xc{0~^) one. This is partially due to a significant contribution 
from lower polarization states A = for Xc{^~^) and A = 0, ±1 for Xc{'^~^) meson. Corresponding 
theoretical uncertainties are discussed. 
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I. INTRODUCTION 



It is well known that the exclusive diffractive Higgs production provides a very conve- 
nient tool for Higgs searches at hadron colliders due to a very clean environment unlike the 
inclusive production A QCD mechanism for the diffractive production of heavy central 
system has been proposed by Kaidalov, Khoze, Martin and Ryskin (Durham group) for 
Higgs production at the LHC (see Refs. (H-IH). Below we will refer to it as the KKMR 
approach. In the framework of this approach the amplitude of the exclusive pp — ?• pXp 
process is considered to be a convolution of the hard subprocess amplitude describing fusion 
of two off-shell gluons producing a heavy system g*g* — )■ X, and the soft hadronic factors 
containing information about emission of the relatively soft gluons from the proton lines 
(see Fig. [1]). In the framework of the fcx-factorisation approach these soft parts are writ- 
ten in terms of so-called off-diagonal unintegrated gluon distributions (UGDFs). The QCD 
factorisation is rigorously justified in the limit of very large factorisation scale being the 
transverse mass of the central system M_|_. 




FIG. 1: The QCD mechanism of diffractive production of the heavy central system X. 



In order to check the underlying production mechanism it is worth to replace the Higgs 
boson by a lighter (but still heavy enough to provide the QCD factorisation) meson which 
is easier to measure. In this respect the exclusive production of heavy quarkonia is under 
special interest from both experimental and theoretical point of view 3]- Verifying the 
KKMR approach against various data on exclusive meson production at high energies is a 
good test of nonperturbative dynamics of parton distributions encoded in UGDFs. 

Recently, the signal from the diffractive Xc(0''', l"*", 2"'") charmonia production in the 
radiative 7/^^-1-7 decay channel has been measured by the CDF Collaboration [sj: 
da / dy\y=Q{pp — )• pp{J /ip + 7)) ~ (0.97 ± 0.26) nb. Assuming the absolute dominance of 
the spin-0 contribution, this result was published by the CDF Collaboration in the form: 



da 
dy 



,(Xc(0+)) 



da 



y=0 



{pp — )■ pp{J/ip + 7)) = (76 ± 14) nb. 



Indeed, in the very forward limit the contributions from Xc(l^, 2"*") vanish due to the Jz = 
selection rule 11, 27]. This is not true, however, for general kinematics 0, 0]- In particular, 
it was shown in Ref. jl] that the axial-vector Xc(l^) production, due to a relatively large 
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branching fraction of its radiative decay, may not be negligible and gives a significant con- 
tribution to the total signal measured by the CDF Collaboration. The same holds also for 
the tensor Xc(2^) meson contribution Recent Durham group investigations 10| support 
these predictions. 

The production of the axial-vector Xc{^~^) nieson is additionally suppressed w.r.t. 
Xci^^j 2"*") in the limit of on-shell fusing gluons (with non-forward protons) due to the 
Landau- Yang theorem Such an extra suppression may, in principle, lead to the dom- 
inance of the Xc(2^) contribution over the Xc{^~^) one in the radiative decay channel jof. 
Off-shell effects play a significant role even for the scalar Xc{0~^) production reducing the 
total cross section by a factor of 2 - 5 depending on UGDFs (6|. The major part of the 
amplitude comes from rather small gluon transverse momenta Q± < 1 GeV. This requires 
a special attention and including all polarisation states Xcll"*", 2+). Our present goal is to 
analyze these issues in more detail in the case of tensor charmonium production at the Teva- 
tron, to study its energy dependence and to compare with corresponding contributions from 
scalar and axial-vector charmonia. 

The paper is organized as follows. Section 2 contains the generalities of the QCD central 
exclusive production mechanism, two different prescriptions for off-diagonal UGDFs are 
introduced and discussed. In Section 3 we derive the hard subprocess amplitude g*g* — )■ 
Xc(2''') in the nonrelativistic QCD formalism and consider its properties. Section 4 contains 
numerical results for total and differential cross sections of Xc{0~^, 1^? 2+) CEP and their 
correspondence to the last CDF data. In Section 5 the summary of main results is given. 



II. DIFFRACTIVE pp ppXc{2^) PRODUCTION AMPLITUDE 

The general kinematics of the central exclusive production (CEP) process pp — )■ pXp 
with X being the colour singlet qq bound state has already been discussed in our previous 
papers on Xc(O^) Q and Xc{^~^) m production. In this section we adopt the same notations 
and consider the matrix element for exclusive Xc(2'^) production and its properties in detail. 

According to the KKMR approach the amplitude of the exclusive double diffractive color 
singlet production pp — )■ ppXcj is 

%,t (lit (lit 

where ti,2 are the momentum transfers along the proton lines, go is the momentum of the 
screening gluon, gi2 are the momenta of fusing gluons, and fgJ{xi,Xi,qQ^,qf-i-,ti) are the 
off-diagonal UGDFs (see Fig. OP). 

Traditional (asymmetric) form of the off-diagonal UGDFs is taken in the limit of very 
small x' <C Xioin analogy to collinear off-diagonal gluon distributions (with factorized 
t-dependence) [ll, l^], i.e. 

/.I = R,f!,'\x^,Qf;;,^^')■F^{t,), 

C = R,ff\x2,Qlf,fi')-F^{h), /.^ = (2.2) 
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with a quasiconstant prefactor Rg which accounts for the single logQ^ skewed effect [14 
and is found to be 1.4 at the Tevatron energy and 1.2 at the LHC energy (for LO PDF), 
Qy2t ~ '^i^(5'o,i' ^1/2 1) ^'^^ effective gluon transverse momenta, as adopted in Ref. P.[ll|. 
F/v(t) is the proton vertex factor, which can be parameterized as F]y{t) = exp{bot) with 
bo = 2 GeV~^ ill, or by the isoscalar nucleon form factor Fi{t) as we have done in Ref. j6[. 
Below we shall refer to Eq. as KMR UGDF^. 

Our results in Ref. [g] showed up a strong sensitivity of the KMRS numerical results 
[ll| on the definition of the effective gluon transverse momenta Qf/2 1 ^'^'^ factorisation 
scales /ii^2- This behavior is explained by the fact that for Xc production the great part of 
the diffractive amplitude (12. ip comes from nonperturbatively small qo,t < 1 GeV. It means 
that the total diffractive process is dominated by very soft screening gluon exchanges with 
no hard scale and extremely small x' <^ xi^2- 

In principle, the factor Rg in Eq. (12. 2 p should be a function of x' and xi or X2- In this 
case the off-diagonal UGDFs do not depend on x' and gg,* {or ql/2t)^ their evolution is 
reduced to diagonal UGDFs evolution corresponding to one "effective" gluon. In general, the 
factor Rg can depend on UGDF and reflects complicated and still not well known dynamics 
in the small-x region. 

In order to test this small-x dynamics and estimate the theoretical uncertainties related 
to introducing one "effective" gluon transverse momentum instead of two ones in Eq. (12. 2p . 



m Refs. |6|, m we have used more generalized symmetrical prescription for the off-diagonal 
UGDFs. Actually, it is possible to calculate the off-diagonal UGDFs in terms of their 
diagonal counterparts as follows^ 



C = (^2,gi„/i^) ■ Fiv(t2) , (2.3) 

where 

I _ I 2 _ 2 2 _ ^^^1. 

^1 — ^25 f^O ~ %,t^ ^ ~ 4 ' 

This form of skewed two-gluon UGDFs (12. 3 p is inspired by the positivity constraints for the 
collinear Generalized Parton Distributions [13], and can be considered as a saturation of 
the Cauchy-Schwarz inequality for the density matrix [isj. It allows us to incorporate the 
actual dependence of the off-diagonal UGDFs on longitudinal momentum fraction of the soft 
screening gluon x' and its transverse momentum q^^ in explicitly symmetric way. As will 
be shown below, these symmetric off-diagonal UGDFs lead to results which are consistent 
with the Tevatron data. 

However, trying to incorporate the actual dependence of UGDFs on (small but neverthe- 
less finite) x' we may encounter a problem. The kinematics of the double diffractive process 
pp — )■ pXp does not give any precise expression for x' in terms of the phase space integration 
variables. From the QCD mechanism under consideration one can only expect the general 



^ In actual calculations we use a more precise phenomenological Martin-Ryskin UGDF introduced in 
Ref. 



13|. We are very thankful to L. Harland-Lang for a discussion on this point. 



^ For diagonal distributions without explicit scale dependences the /Iq, /i^ arguments must be omitted. 
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inequality x' <^ xi,2 and upper bound x' < go,t/ since the only scale appearing in the left 
part of the gluon ladder is the transverse momentum of the soft screening gluon go,t- 

To explore the sensitivity of the final results on the values of x', staying in the framework 
of traditional KKMR approach, one can introduce naively x' = S,- <io,t/ \/s with an auxiliary 
parameter ^ fsj. In our earlier papers we considered the limiting case of maximal 

x' (with ^ = 1). However, it is worth to compare the predictions of the underlying QCD 
mechanism for smaller ^ against the available experimental data in order to estimate typical 
x' values. We will analyze this issue in greater detail in the Results section. 



III. HARD SUBPROCESS g*g* Xc(2+) AMPLITUDE 



Projection of the hard amplitude onto the singlet charmonium bound state V^l'^^ is given 
by an 4-dimentional integral over relative momentum of quark and antiquark q = {ki — k2)/2 



[19|,|20 



VjTAki, k,) = V{qq -> Xcj) . ^2-Uku k2)=2n-J2J2^J (^I^ ' ^ 



i,k Lz,Sz 



X 



X 



^l=i,lM) ■ {L = l,W,S = 1,5,1 J, J,) {3t,3k\l)TT{^':i'^l^Vs=i,s.} 



(3.1) 



■q,ciC2 



-9 



J.C1J.C2 



Qi,t - ki^ 



m 



{qi — kiY — "w-- 



J.C2 J-Cl 



7m 



qi,t - k2,t + m 

{qi ~ ^2)^ ~ 



77. 



Here the function (^l=i,lX^ the momentum space wave function of the charmonium, the 
Clebsch-Gordan coefficient in color space is (3i,3/c|l) = 5"^^ / y/Wl, the trace of t-matrices is 
{t'^^t'^^) = 6'^'^'^^/2, and the projection operator Vs=i,Sz for a small relative momentum q 
has the form 



m] 



V2 



■{ki + m). 



(3.2) 



Since P-wave function $l=i,L3 vanishes at the origin, we may expand the trace in Eq. (13. ip 
in the Taylor series around q = 0, and only the linear terms in q'^ survive. This yields an 
expression proportional to 



(27r 



—I 




;'^(i^.)7^'(o), 



(3.3) 



with the derivative of the P-wave radial wave function at the origin Tl'{0) whose numerical 
value can be found in Ref. ill. The general P-wave result (13.1 p may be further reduced by 
employing the Clebsch-Gordan identity which for the tensor Xcj=2 charmonium states reads 

Taking into account standard definitions of the light-cone vectors n"*" = P2/Ecms, n~ = 
Pi/Ecms and momentum decompositions qi = XiPi + qi^t, q2 = 2:2^2 + q2,t and using the 
gauge invariance property (Gribov's trick) one gets the following projection (for any spin J) 



J, ^lU 



q'2vr 



J, flU 



0, 



Vr'{quq2) 



ntn:VjyMul2) = -^^-^^-^Vj^:^^M,q2). 



(3.4) 
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Since we adopt here the definition of the polarization vectors proportional to gluon transverse 
momenta gi/2,t, then 

K7T(^i.*'^2,t) ^0, gi,t^O, or g2,t^0. (3.5) 

It shows that gluon transverse momenta are necessary to get a nonzeroth diffractive cross 
section. 

Summarizing all ingredients above, we get the vertex factor g*g* — )■ Xc(2^) in the following 
covariant form 



T/C1C2 



7^'(0)e 



(A) 

per 



MttN^ MMl{qiq2) 



(3.6) 



?2,t(?l,t?2,* + llt1l,t) + i^lPl - X2P2)i(ll,t(l2,t - <l2,t<lit) + {<ll,tq2,t){Xlp{ - X2P^2){(ll,t " ^2 

2ql^Xip{ql^ - 2ql^X2P^qlt + 2{qi^tq2,t){xiPiq2,t + X2P2Qi,t) + -^iQi,tq2,t)iPiP2 + P2P1) 
Polarization tensor of Xcj=2 satisfies the following relations (see e.g. Ref. (22| ) 
P%,(A) = P^e^,(A) = 0, e^,(A) = e,^(A), e^^(A) = 0, e^,(A)e^-*(A') = 5xx' 



eM-(A)ep/(A) = ^M^pM,, + \m^Mup - ^M^^Mp^, 



M, 



\=0,±1,±2 

One can check that it may be represented in the following general form 



M2 



;±|A|) 



76, 



^0|A| <?^3 + 



9fj.u 



P P 

M2 



(3.7) 



-2'^2|A| \i[ntn2 + «] ± [n'(n'{ - n'^rf^] 

where ni,2,3 are light-like basis vectors satisfying n'^n'^gpi, = gap (with tiq = P^/M), and 
A = 0, ±1,±2 are the Xc(2^) meson helicities. To our best knowledge, there is no ex- 
plicit decomposition of the meson polarisation tensor e^;y(A) in terms of basis vectors Ui like 
Eq. (13.71) in the literature. In practical calculations below it is convenient to use it in a 
different representation: 



e..(A) = ^(2-|A|)(l-|A|) 



P^Pu 

M2 



+ ^(2-|A|)(l-|A|)n;;n3^ + 



+ ^A(1 - |A|)K< - n^^] + ^^|A|(1 - |A|)[nK + n'2<] + 
+ l\{2 - lADKn;^ + n^] + U\\\{2 - \\\)[ni^n^, + n^n^] . 



Similarly to what has been done for Xc(l^) production in Ref. [8|, in the c.m.s. frame we 
choose the basis with coUinear 113 and P vectors (so, we have P = {E, 0, 0, Pz), Pz = |P| > 0) 
as a simplest one 

1 



n 



/3 



(0,1,0,0), = (0, 0, 1, 0), n 



M 



|P|,0,0,i?), 



VE^ - M2. (3.8) 
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FIG. 2: Coordinate basis in the center- of -mass system of incoming protons pi,2- 



Note, that we choose n2 to be transverse to the c.m.s beam axis (see Fig. [2]), while rii, 
are turned around by the polar angle if) = [0 ... vr] between P and the c.m.s. beam axis. In 
the considered basis {rii, n2, n^} we have the following coordinates of the incoming protons 

Pi = —{I, -sintp, 0, cosip), p2 = —{I, sin^, 0, -cosip) . (3.9) 
The gluon transverse momenta with respect to the c.m.s. beam axis are 

qi,t = (0, Qltcosip, Ql^siniP), g2,t = (0, Qlt cos -Q\, Q^^^sin^), 

where Q^/2t' '^Q\ ^"^^ the components of the gluon transverse momenta in the basis with 
the z-axis collinear to the c.m.s. beam axis. 

From definition (13. 9 p it follows that energy of the meson and polar angle ip are related to 
covariant scalar products in the considered coordinate system as ^ 

E=^l}n±lPlEl^ cosi.= ^lin_(^^ sini. = ^If^A^^. (3.10) 



Furthermore, we also see that from qi = Xipi + qi^t, 92 = 2:^2^2 + Q'2,t and gi + g2 = -P we have 

+ I P I cos '?/' E—\P\cosip 
Xi = j= , X2 = 1= . (3.11) 

Relations (I3.10p and (13. lip show that the interchange of proton momenta pi 4-> p2 is equiv- 
alent to the interchange of the angle ^ ■i/' ± vr, i.e. sinip ^ — sinip and cosip ^ — cosifj 
simultaneously. The last permutation also provides the interchange of the longitudinal com- 
ponents of gluons momenta Xi 0x2. 

Conservation laws provide us with the following relations between components of gluon 
transverse momenta and covariant scalar products 



Qit = ipi ■ / ' Q2,t = ■ , , Qt = psn sign 

IPIsm-?/' ' |P|sm?/^ \Pt\ 

P2 = -|P,|2 = -|p|2sin2 V; = qft + qj^ + 2(gi,tg2,t), Qy2,t = -|qi/2,t|', 



where \Pt\ = \P\\smilj\ is the meson transverse momentum with respect to the 2;- axis. 
The appearance of the factor sign(Q^) guarantees the apphcabihty of fl3.12p for positive 
and negative Q^. Note that under permutations qi^t i?2,t imphed by Bose statistics the 
components interchange as Qff ^ Qf t ^^nd ^ —Qt- Iii our notations the quantity sinip 
plays a role of the noncollinearity of meson in considered coordinates. A straightforward 
calculation leads to the following vertex function in these coordinates 



6M'^i\X\{ql^ - ql^) sign(Qn { I [qi,t x ^2,t] x ni| (1 - |A|) sign(sin V^) sign(cos?/;) + 
2 |[qi,t X q2,i] X nsl (2 - |A|)} - [iq^ql, + (ql, + qlt)iqi,tq2,t)] {3M'{cos^^ + 1)A(1 - |A|) + 
6MEsm{2^) A(2 - |A|) sign(sinV') sign(cos V-) + ^6 (M^ + 2E^) sin^ tp {1 - \\\){2 - \X\) 



where 



[qi,t X q2,t] X nil = \ qhqlt - (^i,t?2,t)^ I cost/'I, 



E 



|[qi,t X q^A X nsl = —^qltqlt - {quq2,tf \ sinV'!- 

The amplitude fl3.12p explicitly obeys the Bose symmetry under the interchange of gluon 
momenta and polarizations due to resulting simultaneous permutations cosif) o —cos?/', 
sin-^ ^ —smip and Ql O —Q\- 

It follows from the conservation laws that 

(lit + p'u = -Qot, q2t + P2t = 9ot, Pt = -{p'lt + P2t) 
Let us consider first the limit of the "coherent" scattering of protons p'u = = Pt, so 

Qit = -{pt + qot), q2t = -{pt - qot), Pt = -2pt, p^ = o. (3.13) 

The production vertex fl3.12p in this limit has the form 



VjTlMt^M = 2^/^"^/ ...' . ...... , ^2. .2 X 



SMnNc M{M^ - 2{j>'l + g2j)2 
l2M'i\\\ ql,ql[{l - |A|) cos^ + ^ (2 - |A|) sinV^} + (3.14) 



[pI + {Qot? - {qltf] {3M2(cos2 i, + 1)A(1 - |A|) + 

6MEsin(2V') A(2 - |A|) sign(sin^) sign(cos^) + ^6 (M^ + 2^^) sin^ ^ (1 - |A|)(2 - |A|) 



We see that in contrast to the axial- vector case considered in Ref. [8[, the diffractive ampli- 
tude of Xc(2^) production does not turn to zero in this "coherent" limit for pt 7^ 0. 
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In the forward limit pt ^ (which is a particular case of the "coherent" one) the ampli- 
tude turns to zero at any meson rapidities y. Indeed, we have Pt ^ and sin?/' — )■ ±0 and 
the amplitude turns into 



VSM,P. - 0) ^ /^-^.y^JH^i^ X (3,15) 

(1 - \\\){t\[{q'„f - (Ar-] - 2 |A| sign(cos'i;,)|,^„,} . 

The imaginary part of this vertex function turns out to be antisymmetric w.r.t. interchanging 
Qot ^ itti whereas its real part is antisymmetric w.r.t. changing the sign of q^^ or q^^ 
component, i.e. 

'^VTJlMt.qlt.Pt ^ 0) = -^V^]:iMlt.<fot.Pt ^ 0) 

5^^;i?Ate, 4t.vt ^ 0) = -5Ry;i|,(-go^„ qi.vt ^ 0) = -mrr=2Mt^ -^It^vt ^ o) . 

Since in this case qu = —qot, qit = qot in the forward limit, then the double integral in the 
diffractive amplitude has an antisymmetric integrand and turns to zero in the symmetric 
limit 

Mp,^o ~ Fi{ti)Fi{t2) / dqQ^dq'ot 6 = ^■ 

(3.16) 

This explicitly confirms the observation made in Refs. js], 23 

Very recently, when our paper was almost complete, a paper by L. Harland-Lang, V. 
Khoze, M. Ryskin and W. Stirling (HKRS) [lo[ appeared where the hard subprocess ampli- 
tudes gg — 7- Xc{J^) (based on formalism by Kuhn et al for 7*7* — Xc{J~^) [24]) including 
the gluon virtualities were listed for different spins including the tensor Xci'^'^)'- 

^7=^' = ^ [3M'{qi,tq2,t) - {qi,tq2,t){ql + qlt) ' KA] ' (3-17) 

VfSx = -^Pi,.P2,ae'"'-'e, [iq2,tUlt - iqi.Ult] , (3.18) 

^HKRs ^ + 2iq,,q,,t)pi,,P2,a\ , (3.19) 
where the constant prefactor is 



1 / 6 

The first amplitude Vf^^^^ ^^TH} is the same as the expression obtained in Ref. |Q] (up 
to a factor of 2 coming from different normalisations of the hard part n^n~Vj^^,^ in our case 
and {2/ s)piP2^,j^fj^i, in Ref. [lo|), where the major role of the gluon virtualities in the hard 
subprocess amplitude of quarkonia production was claimed to be crucial. In particular, it 



We are grateful to V. A. Khoze for helpful discussions of this problem. 
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was shown that an account of the gluon virtuahties reduces the previous KMRS result in 
Ref. [HI] for on-mass-shell gluons Vj^^^^ ~ {qi,t<l2,t) by a factor of 2 - 3. 

The second amplitude, V^j^!^^^, looks different from our previous result, obtained in 
Ref. jsf. However, one can directly check that the difference between the amplitudes f l3.18p 
and (2.12) in Ref. jsj turns to zero when fixing the coordinates in the c.m.s. frame of refer- 
ence as in Eq. fl3.9p (see also Fig. [2]) and the meson polarisation vector e'^ with the basis as 
in Eq. fl3.8p . Due to the covariant structure of these amplitudes, the last observation means 
that they are the same in any frame of reference. The calculations proving this equality are 
rather involved, and we do not show them explicitly here. 

Very similar situation holds for Xc(2^) production amplitudes. Namely, the amphtudes 
f l3.19p and (13. 6p turned out to be the same under fixing the coordinates as in the previous 
section. Therefore, under the kinematical relations our results for the hard subprocess 
amplitudes are in complete agreement with the corresponding HKRS results. Let us now 
turn to the discussion of numerical results. 



IV. NUMERICAL RESULTS 

Results for the differential cross sections da/dy{y = 0) of the diffractive Xc(0^,l^,2+) 
meson production at the Tevatron energy W = 1960 GeV for different UGDFs are shown 
in Table [H In the last column we show the results for the expected signal in the J/tp + j 
channel summed over all Xc spin states (and all polarisation states of Xc(l^; 2+) mesons) 



da. 



obs 



dy 



, ^ E K^2o{si^) jmxc{.n ^ + 7)^^ , (4.1) 

y=o dy y=o 



which can be compared with the corresponding value measured by the CDF Collaboration 
Q: da^''P/dyL=Qim) pp{J/ip + 7)) ~ (0.97 ± 0.26) nb. 

In Refs. |lOl . Ill| is was assumed that the NLO corrections factor -f^NLO in the g*g* — x 
vertex is the same as in the x ^ 99 width implying that |Vjp ~ r{x 99). In general, 
such corrections depend on spin of qq resonance. So, the diffractive cross section for each 
Xc{J^) has to be multiplied by not necessarily the same factor K^^q, as shown in Eq. (14. ip ^. 
This can be done, however, only for 0"*" and 2^ states, and the corresponding NLO QCD 
radiative corrections are well-known (25| : 

iriSo = 1 + 8.77^^ ^ 1.68, i^iSo = 1 - 4.827^^ ^ 0.63 . (4.2) 

71 71 

Due to the Landau- Yang theorem the decay of the axial vector charmonium l"*"^ to on-shell 
gluons is forbidden, and there are no reliable calculations of the NLO QCD corrections to 
its coupling with off-shell gluons. In the following we take naively -R^nlo ~ This leads to 
an additional uncertainty of the model predictions. 

As has been claimed in Refs. 0, the absorptive corrections are quite sensitive to the 
meson spin-parity. This was studied before in the context of scalar and pseudoscalar Higgs 



See Ref. [l^l for discussion of extra uncertainties coming from NNLO and higher order corrections. 
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production in Ref. [2[. We adopt here the following effective gap survival factors, calculated 
in Ref. [lo[ for different spins including eikonal and so-called enhanced contributions: 

(5cff(Xc(0+))) ^ 0.033, (^eff(Xc(l+))) ^ 0.050, (^,^ff(Xc(2+))) ^ 0.073 . (4.3) 

The contribution of the scalar Xc(O^) CEP, which was initially assumed to be the dom- 



inant one is reduced by a very small branching ratio of its observable radiative decay 
js], In turn, the strong suppression of the Xc(l^) central production in both the on-mass- 
shell limit of fusing gluons (due to Landau- Yang theorem 26|) and the forward scattering 
limit of outgoing protons (due to the so-called = selection rule [III |2^) may be par- 
tially compensated by its much higher branching ratio to the observed J/?/' + 7 final state |8j]- 
Analogously to the axial-vector case, the suppression of the tensor Xc{'^~^) CEP is likely to 
be eliminated by its large decay branching ratio 10(|, and the resulting value of the radiative 
decay signal is under our special interest. 



TABLE I: Differential cross section da^^/dy{y = 0) (in nb) of the exclusive diffractive produc- 
tion of Xc(0"'", 1^, 2+) mesons and their partial and total signal in radiative J/ifj + 7 decay channel 
daji^^/dy{y = 0) at Tevatron for different UGDFs, cuts on the transverse momentum of the glu- 
ons in the loop {qo,t) and different values of the auxiliary parameter ^ controlling the characteristic 
x' values in the symmetric skewed UGDFs prescription ()2.3I) (denoted as "sqrt"). NLO skewedness 
factor R^^'^ = 1.3 for the KMR asymmetric prescription ()2.2p (denoted as "Rg"), NLO correction 
factors ()4.2p and absorptive correction factors ()4.3p are included. Contributions from all polarisations 
are incorporated. 



skewed UGDF 
prescription 




Xc(0+) 


Xc(l+) 


Xc(2+) 


ratio 


signal 

ddoba 

dy 


dy 


<i'^j/ip-f 
dy 


dy 


do'.l/ip-f 
dy 


dy 


d(y jj^-^ 
dy 


1+/0+ 


2+/0+ 


GBW f30j, "sqrt" 


1.0 


13.2 


0.15 


0.01 


0.003 


1.96 


0.38 


0.02 


2.5 


0.5 


0.3 


12.8 


0.15 


0.04 


0.01 


1.39 


0.27 


0.07 


1.8 


0.4 


lin KS [31], 




32.6 


0.37 


0.20 


0.07 


0.53 


0.10 


0.2 


0.3 


0.5 


lin KS [31], "sqrt" 


1.0 


17.2 


0.19 


0.12 


0.04 


0.61 


0.12 


0.2 


0.6 


0.4 


nlin KS [3JJ, "sqrt" 


1.0 


12.6 


0.14 


0.07 


0.02 


0.39 


0.08 


0.1 


0.6 


0.3 


0.3 


20.6 


0.23 


0.10 


0.03 


0.57 


0.11 


0.1 


0.5 


0.4 


0.05 


36.0 


0.41 


0.13 


0.04 


0.84 


0.16 


0.1 


0.4 


0.6 


KMR [13], GRV94HO, 
(gg«/)2 = 0.72 GeV^ 




13.5 


0.16 


0.06 


0.02 


0.23 


0.05 


0.1 


0.3 


0.2 


KMR [13], GRV94HO, 
(gg«t)2 = 0.36 GeV^ 




37.9 


0.43 


0.14 


0.05 


0.75 


0.15 


0.1 


0.3 


0.6 


HKRS result [10] 
(gg"/)2 = 0.72 GeV^ 




27.1 


0.31 


0.72 


0.25 


0.95 


0.19 


0.8 


0.6 


0.7 



As it was discussed in Ref. [6|, the dominant contribution to the diffractive CEP of Xc(0''') 
comes from nonperturbative values of the gluon transverse momenta < 1 GeV. In order to 
estimate the role of small qt in the central production of Xc(1^5 2+) and related theoretical 
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uncertainties we use different UGDFs known from the literature (for details see Refs. BHil). 
Among them there are perturbatively modeled KMR UGDF [l|, |l2|, |l3|, which include the 



Sudakov form factor, as well as GBW ^] and linear /nonlinear Kutak-Stasto (KS) j3l| 
UGDF models which by construction can be used for any values of the gluon transverse 
momenta. 

In the last row of Table^we show the HKRS results for partial cross sections extracted 
from their original paper [lOj]. These cross sections were calculated by the HKRS at some 
small energy scale and subsequently extrapolated up to the Tevatron energy assuming a 
Regge type energy dependence. 

By direct calculation at the Tevatron energy with the same KMR UGDFs, but with- 
out imposing any arguments beyond the QCD framework (like Regge scaling, for ex- 
ample), we get the observable J/ip'~i cross section daj/^^/dy{y = 0) ~ 0.6 nb, which 
is close to HKRS result ~ 0.7 nb, but somewhat lower than the present CDF result 
da^^P / dy\y=Q{pp — )■ pp^J/ip + 7)) ~ (0.97 ± 0.26) nb. However, this demanded to incorpo- 
rate physics below HKRS cut-off on gluon transverse momentum in the loop integral (12. ip 
i^of)'^ = 0.72 GeV underlining the importance of nonperturbative contributions of small go,« 
in the QCD mechanism under consideration. Relations between different Xc's obtained in 
Ref. are not reproduced as well. Our result is highlighted in bold in Table [H 

The reason of such a discrepancy in both the normalization and relative contributions of 
different spins to the observable rate is the on-shell approximation 3> 1 hard 
9*9* Xc{J~^) adopted by HKRS^. As was firstly noticed in Ref. |6?], the gluon off-shellness 
plays an important role in the central exclusive Xc{0~^) leading to a strong reduction of the 
cross section depending on UGDF model. Now, we observe that the same effect significantly 
affects the ratios between different spin contributions. 

Note that at zeroth meson rapidity y = a. significant part of the cross section comes from 
lower polarisation states in the center-of-mass frame A = {xc{^~^)) and A = 0, ±1 (xc(2^)). 
In the total (integrated over y) cross section the maximal helicity contributions, however, 
strongly dominate. We leave a more detailed investigation of the polarisation effects for a 
separate publication. 

Relative contributions of Xc(0^, 1^,2+) CEP to observable signal {J/"^ + 7) require an 
additional discussion. Last PDG updated set of branching ratios for charmonia radiative 



decays is [28 



BR(xc(0+) ^ J/iJ + 7) = 0.0114, 
BR(xc(l+)^ J/V^ + 7) = 0.341, 
BR(xc(2+) ^ J/tp + 7) = 0.194. 

Furthermore, as one can see in Table HI despite of larger branching ratio in the axial- vector 
case the observable signal from the Xc{^~^) CEP occurs to be smaller than that from xd'^'^) 
for UGDFs enhanced at sufficiently small nonperturbative qt (in particular, for the Kutak- 
Stasto (KS) and GBW UGDFs) due to an additional suppression of the g*g* — )■ Xc(l^) 
subprocess vertex at small qi/2,t- For the GBW UGDF the Xc{^~^) contribution is strongly 
suppressed whereas the Xc(2^) contribution turned out to be larger than the Xc(O^) one. All 



^ We are most thankful to L. Harland-Lang for the helpful discussions and correspondence on the topic. The 
exchange of the Fortran codes between us helped a lot in cross-checks of our calculations and understanding 
the various sources of discrepancies. 



12 



do/dy(y=0), nb 









X,(0') 




-- x,(l1 




----X,(2*) 


— \ 

1 ' 1 ' 1 ' 1 ' 1 ' 1 


x 

' 1 1 1 




0.2 0.3 0.4 0.5 0.6 0.7 


0.8 


1 ' 1 
0.9 1.0 



Q', GeV' 

t.cut 

FIG. 3: Dependence of the differential cross section da^^/dy{y = 0) 0/ Xc(0"'", 1^, 2+) CEP on 
the infrared cut-off on small effective gluon transverse momentum Qf^^ for the KMR UGDF with 
GRV94HO (Rg = I.3J. Absorption effects are not included here. Arrow points to the HKRS cut-off 
0.72 GeV^ 0y. 

the UGDF models under considerations lead to somewhat underestimated observable signal 
at Tevatron < 0.6 nb, however, it can be still reliable within relatively large theoretical 
uncertainties of the gap survival factors and QCD mechanism under discussion. In the case 
of the KS model, contributions from XciX^) state are found to be stronger suppressed than 
in KMR model. Measurements of separate spin contributions thus would help to distinguish 
between different UGDF models. 

TABLE II: Integrated over full phase space (bare) cross sections (in nb) for the central exclusive 
Xc(0+,l+,2+) production at RHIC, Tevatron and LHC energies. Infrared cut-off for the KMR UGDF 
(in "i?^" -prescription) is taken to be (Qj"*)^ = 0.72 GeV^. We take Rg to be equal 1.3 at all three 
energies. Absorption effects are not included here. 



Xc 


UGDF 


RHIC 


Tevatron 


LHC 


Xc(0+) 


nlin. KS, i = 0.3 


108 


3569 


23260 


KMR, GRV94HO 


43 


2270 


30720 


Xc(l+) 


nlin. KS, i = 0.3 


0.4 


12 


65 


KMR, GRV94HO 


0.2 


7 


87 


Xc(2+) 


nlin. KS, i = 0.3 


2 


44 


209 


KMR, GRV94HO 


0.4 


18 


195 



In the case of the KMR UGDF, we observe quite substantial dependence of the predicted 
observable signal w.r.t. variations of the infrared cut-off on small transverse momenta of 
the gluons in the most internal loop (see Fig. 12]) • From Table H] we see that the shift of Q^"* 
from the value 0.72 GeV^ used in Ref. [13] down to the minimal perturbative scale of the 
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integrated GRV94HO distributions 0.36 GeV^ Q leads to increase of the cross section by 
a factor of about 3, approaching the CDF data. For comparison, decrease of the Qcut from 
1 GeV^ down to 0.36 GeV^ leads to increase of the cross section by a factor of 6. Since we can 
not estimate the nonperturbative contribution coming from below 0.36 GeV^, this allows us 
to conclude that perturbatively motivated KMR UGDF leads to infrared unstable result in 
the case of relatively light charmonium CEP. It is clear that the essential part of the QCD 
dynamics comes from the nonperturbative region of transverse momenta below the HKRS 
cut-off 0.72 GeV^ [lo|. KS and GBW UGDFs allow to incorporate some unknown physics 
even below the minimal GRV scale Qq = 0.36 GeV^, avoiding ambiguities in defining the 
effective gluon momenta. 

Applying the KMR's asymmetrical off-diagonal UGDF according to Eq. (12.21) ( "-Rg" pre- 
scription) in the case of the GBW models we get strongly overestimated observable signal at 
Tevatron, which means that in this case it is crucial to take into account the x'-dependence 
of off-diagonal UGDFs when going deeply into the infrared region of small g^'s. The x'- 
dependent "sqrt" prescription, introduced in Eq. (12. 3p . leads to observable signal, which is 
much closer to the experimental data. 

The "sqrt" prescription, introduced in Eq. (12. 3p . provides an agreement with the data 
(within a factor of 2 in overall theoretical uncertainty between different UGDFs) with the KS 
(with rather small ^ ~ 0.05) and GBW models giving the cross section daobs/dy{y = 0) ~ 
0.5 — 0.6 (see Table[T]). This practically means that the smaller go,t comes into the game, the 
smaller x' w.r.t. q^t/s is required to get the data description, providing one more argument 
about importance of nonperturbative effects in charmonia CEP. The relative contributions 
of different charmonium states in the J/ip + j channel (including absorption effects) in the 
case of, e.g., KS model are found to be: 

Jlin\ .IJNn\ .IJlti\ =1.0.1:0.4. (4.4) 
dy I \ dy I \ dy J 

" / KS \ " / KS \ " / KS 

They are not affected by smaller x' or nonlinear effects in this model. As the normalization 
point we took the contribution of the Xc(O^) meson CEP as was done in Ref. M. 



In Table |T] we also presented results with the linear Kutak-Stasto model based on the 
unified BFKL-DGLAP framework and the nonlinear one based on the Balitsky-Kovchegov 
equation jill]. It turned out that incorporation of the nonlinear effects responsible for the 
gluon recombination in this model reduces the Xc{J^) CEP cross sections by 30-50 %. We see 
that the nonlinear effects play a crucial role in diffractive quarkonia production effectively 
decreasing the characteristic values of x' (controlled by ^). However, reliable predictions 
including the nonlinear effects require the exact knowledge of the triple Pomeron vertex at 
NLLx accuracy, which is yet unknown. 

It is also interesting to compare the diffractive production of Xc states at different energies. 
As an example, in Table [TTl we present the integrated (over full phase space) cross sections 
of Xc{0~^, l"*", 2"*") production at RHIC, Tevatron and LHC energies. The results show similar 
energy behavior of the diffractive cross section for different UGDFs as well as for different 
Xc states. 

Finally, let us turn to differential distributions. In Fig. H] we show the differential cross 
section da/dy in rapidity y for all Xc states. In this figure and in the following, all helicity 
contributions for Xc(l^; 2^) CEP are taken into account. Here and below we show only bare 
CEP cross sections for GBW, KS and KMR UGDFs. In the last case, we present the results 



computed with the HKRS cut-off parameter 0.72 GeV^ [10[. We see that the shape of the 
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dc/dy, nb 



dc/dy, nb 

20- 




FIG. 4: Distributions da^^ / dy in rapidity of xd^^) (left panel), Xc(l^) (middle panel) and Xc{'^^) 
(right panel) mesons for different UGDFs at the Tevatron energy y/s = 1.96 TeV. The dash-dotted 
line corresponds to the KS UGDF fsi l in the symmetrical "sqrt" -prescription with ^ = 0.3, solid 
line - KMR UGDF 0/ with Rg = 1.3, (Q^*)^ = 0-72 GeV^ and GRV94HO PDF [3^], and short- 
dashed line represents result with the GBW UGDF fSui I (^ = 0.3). Absorption effects are not 
included here. 



dn/dt, nb 



dn/dt, nb 




0.0 0.2 



FIG. 5: Distribution in ti^2 o/Xc(0^) (left panel), Xc(l^) (middle panel) and XcC^^) (right panel) 
for meson CEP for different UGDFs. The meaning of curves here is the same as in Fig. [7} 



curves is rather similar, however, they have substantially different maxima. The biggest 
cross section for the Xci^^, 2"^) states is obtained with the KS UGDF, whereas for Xc{^~^) 
the KS and KMR UGDFs give quite similar cross sections. 

In Fig. 13 we present corresponding distributions in t = ti or t = ^2 (identical), again for 
different UGDFs. Except of normalisation the shapes are rather similar. This is because of 
the ti and t2 dependencies of form factors (describing the off-diagonal effect) are taken the 
same for different UGDFs. 

In Fig. [6] we show the correlation function da/d^ in relative azimuthal angle $ between 
outgoing protons for different Xc states. The shapes of the distributions are somewhat 
dependent on UGDFs. It is interesting to note here that the KS and KMR UGDFs lead 
to very similar angular dependence of da/d^ for all Xc states. In the case when energy 
resolution is not enough to separate contributions from different states of Xc (Xc(O^), Xc{^^), 
Xc{'^~^)), which seems to be the case for Tevatron, the distribution in relative azimuthal angle 
may, at least in principle, be helpful. 

The fact that the angular distributions are not simple functions (like sin $, cos $) of the 
relative azimuthal angle between outgoing nucleons is due to the loop integral in Eq. (12. ip 
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dcj/do, nb 

2000- 



da/cHi, nb 



dn/do, nb 

30- 



30 60 90 120 150 180 30 60 90 120 150 180 30 60 90 120 150 180 

$ (J) o 

FIG. 6: Distribution in relative azimuthal angle $ between outgoing protons for Xc{^^) (left panel), 
Xc(l^) (middle panel) and XcC^^) (right panel) meson CEP for different UGDFs. The meaning of 
curves here is the same as in Fig. [7} 



which destroys the dependence one would obtain with single fusion of well defined (spin, 
parity) objects (mesons or reggeons) 



V. CONCLUSIONS AND DISCUSSION 



Our results can be summarized as follows: 

We have derived the QCD amplitude for central exclusive production of tensor Xci'^'^) 
meson. This amplitude vanishes in the forward limit of outgoing protons, as demanded 
by the = selection rule. Our numerical results show the importance of non-forward 
corrections, including all polarisation states of Xc(2^) and nonperturbative contributions 
to the Xc{'^~^) CEP. Inclusion of all the ingredients leads to a noticeable contribution of 
the Xc{'2~^) meson in the observable radiative decay channel depending on UGDF. We have 
observed the importance of the A = state Xc{^~^) CEP and A = 0, ±1 states for Xc{'2~^) 
CEP at ?/ ~ 0, whereas the total CEP cross section is dominated by maximal helicity 
contributions. 

The main contribution to diffractive charmonium production comes from small gluon 
transverse momenta Qt < 1 GeV leading to quite substantial sensitivity of the corresponding 
cross section on the infrared cut-off in perturbatively modeled KMR UGDF. Alternatively 
one could use UGDFs like Kutak-Stasto and GBW models, which by construction can be 
used for any values of the gluon transverse momenta. 

We have tested the symmetrical prescription for off-diagonal UGDFs, following from pos- 
itivity constraints and incorporating x, qt dependence of both participating gluons, against 
the present CDF experimental data. A rather good quantitative agreement with the CDF 
data on charmonium CEP in the radiative decay channel is achieved with the nonlinear 
Kutak-Stasto UGDF model giving the cross section daobs{,J I^Jl) / dyiy = 0) ~ 0.6 nb with- 
out imposing extra normalisation conditions beyond the QCD framework. Such a description 
is achieved by incorporating very soft screening gluons with x' ~ 0.1 ■ go,t/ a/s- We have also 
calculated total cross sections of Xc CEP at different energies (RHIC, Tevatron and LHC), 
as well as differential distributions in three phase space variables t, 

Overall theoretical uncertainty of the QCD mechanism under consideration is rather 
high but hard to estimate due to large unknown nonperturbative contributions coming 
into the game and not well known higher-order QCD corrections to the hard subprocess 
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g*g* — Xc (especially, in the axial- vector case). Also, absorptive corrections may depend on 
UGDF used in the calculation, and there is no reliable estimation of such a sensitivity in 
literature. In the present paper we kept the strategy to study different distinct options and 
analyze the sensitivity of the final results with respect to the UGDFs choice, prescriptions 
for skewed UGDFs, nonperturbative cut-off parameter and characteristic x' variations, etc. 
Then a comparison with experimental data would allow to select the most reliable option. 
However, we observe a variety of such "good" options, namely, description of the data (with, 
however, pretty large theoretical uncertainties related, in particular, with unknown NLO 
corrections) can be, in principle, achieved for a few UGDFs (GBW, KS and KMR UGDFs, 
see Table [T]). Each of them pick up some essential QCD dynamics. Further constraints can, 
in principle, be settled by experimental measurements of separate Xc{J~^) contributions, the 
energy dependence of the cross section and the shapes of differential distributions. 
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